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Abstract. In recent experiments the irradiation of insulators of perovskite type 
with swift heavy ions under glancing incidence has been shown to provide a unique 
means to generate periodically arranged nanodots at the surface. The physical origin 
of these patterns has been suggested to stem from a highly anisotropic electron 
density distribution within the bulk. In order to show the relevance of the electron 
density distribution of the target we present a model calculation for the system Xe +23 
— > SrTi03 that is known to produce the aforementioned surface modifications. On 
the basis of the Lindhard model of electronic stopping, we employ highly-resolved ab- 
initio electron density data to describe the conversion of kinetic energy into excitation 
energy along the ion track. The primary particle dynamics are obtained via integration 
of the Newtonian equations of motion that are governed by a space- and time- 
dependent friction force originating from Lindhard stopping. The analysis of the local 
electronic stopping power along the ion track reveals a pronounced periodic structure. 
The periodicity length strongly varies with the particular choice of the polar angle 
of incidence and is directly correlated to the experimentally observed formation of 
periodic nanodots at insulator surfaces. 



PACS numbers: 61.80.Jh,61.80.Az,61.82.Ms 



Keywords: electronic stopping, swift heavy ions, incidence, perovskite, insulator, thermal 
spike Submitted to: J. Phys. : Condens. Matter 



Ab-initio calculation of electronic stopping power along. 



2 



1. Introduction 

The irradiation of a solid surface with a ~100 MeV ion beam leads to the formation of 
cylindrical tracks within the target material (see e.g. [HISIE]). This phenomenon can 
be explained in terms of kinetic excitation of electrons along the path of the penetrating 
ion. These electronic excitations are directly connected with structural changes within 
the damage zone which may also extend to the very surface, for instance, manifesting 
themselves in the formation of nanodots, for instance. 

In the case of normal incidence of the ion beam the nanodots are statistically distributed 
over the entire irradiated area [U El El E] due to the random distribution of the ion 
impact sites. In recent experiments [8], the irradiation of perovskite crystals (ABO3) 
under glancing incidence has been shown to provide a unique means to generate chain- 
like, periodic nanodots at the surface, the extension of which varies in the range of a 
few hundred nanometers up to some micrometers depending on the particular angle of 
incidence. This finding has been ascribed to the specific anisotropic electron density of 
the perovskite structure. The nanodots can be imaged by atomic force microscopy [9], 
for instance, and may be regarded as a fingerprint of the electron density projected onto 
the surface. 

Although the detailed microscopic mechanisms of this kind of defect formation are not 
entirely understood, the anisotropy of the electron density inducing a discontinuous 
kinetic energy dissipation of the ion appears to play a key role. For many practical 
purposes the energy loss per track length can be calculated using SRIM pHl E] or 
equivalent standard codes. However, these codes are not capable to account for the 
detailed electronic structure and are limited to perpendicular incidence. In addition, 
these codes do not provide any space-resolved information on energy losses inside the 
crystal, which are supposed to be essential in our case. 

Therefore this paper aims at a detailed investigation of the spatial electronic energy loss 
under ~100 MeV glancing incidence irradiation. For the model system Xe +23 — > SrTiC>3 
which has been shown to produce the aforementioned periodic nano patterns, we 
combine ab-initio electron density computations with trajectory calculations in order 
to study the electronic stopping dynamics along the ion track for different polar and 
azimuthal angles of incidence. 

2. Model 

For the irradiation of SrTi0 3 (100) with Xe +23 ions of ~100 MeV of kinetic energy 
electronic stopping dominates over nuclear stopping by two orders of magnitude. In this 
energy regime electronic stopping occurs via direct inelastic projectile-electron collisions 
and is usually treated in terms of the Lindhard model [12j yielding an electronic stopping 
cross section S e = (3 ■ v p with a material parameter (3 depending on the specific target- 
projectile combination and v p denoting the projectile velocity. For a target containing 
different elements i with individual /3^'s the effective electronic stopping cross section 
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Figure 1. [figure and caption reprinted from Ref. [8] ] DFT calculation of the electron 
density of SrTiC>3: The grey shading represents electron density. Atom positions (apart 
from the central titanium atom) are also visualized. The arrow indicates a possible 
projectile trajectory defined by the azimuthal angle ip as well as the incidence angle 0. 



S e ,es is given by 

S e ,« = P«-V p = Y l U>®P®-V p (1) 

i 

with relative weight factors mirroring the stoichiometry of the target. The electronic 
energy loss per track length (stopping power) of the projectile along its path through 

the solid then reads 
dE 

-^\e = N e (f p (t))-^-\m\ , (2) 

where N e (r p (t)) denotes the electron density of the target at the time- dependent posi- 
tion f p (t) of the penetrating ion. It should be mentioned that for the system studied 
here the initial charge of q = +23 is close to the equilibrium charge state during the 
passage of the ion through the target [131 E] • 

The electron density N e (r) constituting the first essential input parameter in Eq. (j2J) is 
calculated for one unit cell of the perovskite using the Abinit package [15J with pseu- 
dopotentials from the Fritz-Haber-Institut [16]. Employing a cutoff energy of 96 Hartree 
with a sampling of 512 fc-points on a cubic grid, the calculations have been carried out 
within the GGA approximation using the PBE functional [17] . The lattice constant of 
SrTiC>3 has been determined as 3.905 A by minimization of the energy. It should be 
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Figure 2. Average electronic stopping power (dE/dx) e i for different irradiation 
energies (10 MeV, 50 MeV, 90 MeV) and a polar angle of 0.5 °. The step width 
Atp has been set to 0.1 °. 



noted that ab-initio calculations using pseudopotentials only give the electron density 
of the valence electrons. The core electrons, however, are strongly localized around the 
nuclei which limits their influence to a very small fraction of the volume. In the energy 
range we are considering here, the projectile comes close to the nuclei only very rarely. 
Therefore consideration of the core electrons would lead to minor quantitative changes 
of our result, only, and no qualitative changes at all. 

The electron density N e (r) resulting from this ab-initio computation is shown in 
Fig. [TJ The majority of the electrons are localized around the oxygen atoms. Moreover, 
the electron density in the TiC>2 plane clearly exceeds the one observed in the SrO plane. 
Thus, it becomes directly obvious that N e (r) varies on a sub- Angstrom scale and hence 
should not be treated as a homogeneous electron gas. 

The second parameter to be inserted in Eq. ([2]) is the trajectory r p (t) and its first 
time derivative. The trajectory is calculated by numerical integration of the Newtonian 
equations of motion for the projectile. In view of the fact that energy dissipation is 
governed by electronic stopping, in a first-order approximation no interatomic potentials 
have to be incorporated into the equations of motion. Thus, the only interaction between 
the primary ion and the target is an effective friction force according to Eq. (T5]) that 
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Figure 3. Stopping power as a function of ion energy for different polar angles of 
incidence. For each ion energy the stopping power has been averaged over the entire 
range of the azimuthal angle. The dashed line represents the calculated energy loss 
using SRIM [10]. 

originates from Lindhard stopping. During the time integration of the ion track, the 
energy loss per numerical time step is output for further analysis. 

3. Results 

The computational procedure as outlined in the previous section has been performed 
for the irradiation of SrTi03 with Xe +23 ions in the energy range from (10-100) MeV. 
For each bombarding energy the azimuthal angle of impact has been varied. The crystal 
symmetry of the target allows to limit the azimuthal variation from to 45 degrees. 
The impact point is chosen as illustrated in Fig. [1] and kept constant for all calculations 
presented here. However, preliminary computations have shown a negligible influence 
of the choice of impact point particularly with regard to the results presented here. 
Figure [2] shows the dependence of the electronic stopping power averaged along the 
trajectory on the azimuthal angle of incidence for exemplary irradiation energies of 
10, 50 and 90 MeV. The polar angle has been chosen as = 0.5 ° with respect to 
the surface plane as indicated in Fig. [TJ Considering the obtained stopping powers we 
observe rather moderate fluctuations (a ~ 3%) around the mean values of 21.1 keV/nm 
(90 MeV), 15.3 keV/nm (50 MeV) and 6.7 keV/nm (10 MeV), respectively. The obtained 
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finestructures in the angular dependence of the stopping power are directly related to the 
inhomogeneous electron density distribution and exhibit identical curve characteristics 
independent of the ion energy. The pronounced peaks at 19°, 27° and 34° correspond to 
low-indexed crystallographic directions for the chosen impact point. Close to the [001] 
direction minima in the electronic stopping are due to the presence of the Sr atom. The 
ratio of the amplitudes between the different projectile energies scales with \J~E due to 
the employed Lindhard stopping cross section. 

In order to verify the accuracy of our approach we calculate the stopping power averaged 
over the azimuthal degree of freedom for different ion energies and compare these values 
with SRIM calculations. The results are depicted in Fig. [3] for various polar angles of 
incidence 9. It becomes directly obvious that the stopping power is largely independent 
of the angular degrees of freedom. In the energy range from (30-80) MeV the computed 
stopping powers are in accordance with SRIM data. For impact energies below 30 MeV, 
however, we observe a systematic overestimation compared to the SRIM results. This 
discrepancy clearly stems from the fact that in this energy regime the electronic stopping 
power is not acurrately described by a V^-dependence as assumed in our case. 
Summarizing up to here, the averaged electronic stopping power appears to be 
independent of the particular choice of angular impact parameters. However, 
experiments cleary reveal the explicit dependence of the geometry of nano surface 
dots (e.g. total chain length, periodicity length) on the particular choice of polar 
angle. Therefore, it is necessary to take a space-resolved look on the energy transfer 
characteristics along the ion track for different polar angles. 

Figure H] illustrates the electronic stopping power along the 93 MeV Xe 23+ ion track 
in SrTi0 3 for 6 = 0.5°, = 1°, = 2° and = 4°, respectively. For the purpose 
of clarity the depicted trajectory length is limited to 100 nm. The azimuthal angle is 
chosen as <p = 7° representing a highly indexed direction. For = 0.5° (see Fig. H] (a)) 
we notice two pronounced peaks at 20 nm and 70 nm track length. In the middle of 
those peaks we detect more thinned-out secondary peaks of approximately the same 
half width contributing about 35 % to the total excitation energy. Each of the peaks 
exhibits a finestructure of 5-shaped peaks naturally resulting from fluctuating electron 
densities along the track. 

At this point we remark that a higher resolved representation of the data would mirror 
the energy deposition during the passage of each unit cell. However, in our case the 
main point of interest is the observed distance between two consecutive areas of increased 
electronic stopping. The periodicity length derived for = 0.5° is approximately 50 nm. 
For a slightly increased polar angle of 1° (see Fig. |4] (b)) the same kind of peaks can be 
identified, however, each kind with a reduced periodicty of approximately 25 nm. For 
even larger polar angles of = 2° and = 4°, respectively, the two types of peaks 
begin to merge leading to a nearly constant enveloping function. Thus, the pronounced 
peak structure as seen before is increasingly smeared out and the periodicity tends to 
disappear. 

These results qualitatively explain the experimental findings that the periodic nanodot 
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Figure 4. Electronic stopping power as a function of track length, i.e. travelled 
distance of the ion from the point of impact. Data has been calculated for a Xe +23 
ion of 93 MeV initial kinetic energy penetrating a SrTi03 target under polar angles 
9 = 0.5° (a), 9 = 1° (b), 6 = 2° (c) and 6 = 4° (d), respectively. The arrows 
represent the periodicity lengths of subsequent peaks of same type as explained in 
detail in the text. 



formation is limited to glancing angles of incidence up to a few degree. 
We are well aware that a more quantitative treatment of defect formation using our 
model must incorporate a realistic description of the transport of the excitation energy 
as well as the coupling of the electronic system to the lattice. This may be achieved 
in terms of a two temperature approach [18J as employed in the thermal spike model 
[2], for instance. However, in our case the anisotropy of the electron density and the 
necessity to properly implement the surface leads to a break-down of the radial symmetry 
usually presumed in standard thermal spike calculations and, thus, to a drastical increase 
in computational effort to numerically solve the corresponding transport equations. 
Calculations of that type are currently under way in our lab. 
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4. Conclusion 

We have presented a model that uses highly-resolved ab-initio electron density data in 
combination with molecular dynamics in order to obtain space-resolved electronic ex- 
citation characteristics along the track of swift heavy ions in matter. The model has 
been applied to the irradiation of SrTi0 3 with swift Xe +23 ions in glancing incidence 
geometry for different ion energies and angular configurations of incidence. 
Results for the average electronic stopping power are in good agreement with SRIM 
and, moreover, do not show a significant influence of the particular choice of the angu- 
lar degree of freedom. However, the specific choice of the polar angle turns out to be 
crucial for the spatial distribution of electronic excitation along the ion track. A detailed 
analysis of the local behaviour of the electronic stopping power along the trajectory of 
the primary particle reveals a pronounced periodic structure that qualitatively explains 
the physical origin of periodic nanodot formation at insulator surfaces. 
Finally, we would like to point out that our approach is applicable to any insulating 
material for which the electron density data can be calculated. 
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